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Transverse thermal depinning and nonlinear sliding friction of an adsorbed monolayer 
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We study the response of an adsorbed monolayer under a driving force as a model of sliding friction 
phenomena between two crystalline surfaces with a boundary lubrication layer. Using Langevin- 
dynamics simulation, we determine the nonlinear response in the direction transverse to a high 
symmetry direction along which the layer is already sliding. We find that below a finite transition 
temperature, there exist a critical depinning force and hysteresis effects in the transverse response 
in the dynamical state when the adlayer is sliding smoothly along the longitudinal direction. 
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Driven lattice systems interacting with pinning poten- 
tials have attracted growing interest recently due to the 
variety of possible nonequilibrium dynamic phases. This 
problem is relevant for many systems including adsorbed 
layers in tribology , charge-density waves [|j and vor- 
tex lattices JTo|-JT2[ in type-II superconductors. In partic- 
ular, in the boundary lubrication problem for interfaces 
1^, the macroscopic sliding friction is controlled by the 
response of the tightly bound lubricant layer between the 
surfaces to the external driving force. It has been shown 
that there are many dynamical phases, and the transi- 
tions between these phases can lead to such macroscopic 
behavior as the stick and slip motion. The essence of the 
physics in the boundary lubrication has been reasonably 
well understood with the study of a zeroth order model 
consisting in an overlayer of interacting atoms adsorbed 
in a periodic potential subject to an external driving force 
However, an interesting aspect of sliding friction 
which has not been addressed so far is how the nonlin- 
ear response of the overlayer is affected by the symmetry 
properties of the underlying periodic pinning. Initially, 
the overlayer can be driven to slide along a high sym- 
metry direction of the surface. In this sliding state, the 
response to an additional force in the transverse direc- 
tion is of interest. It should provide important informa- 
tion on the nature of the sliding state. In addition, the 
transverse force can also generate novel nonequilibrium 
dynamic phases of the overlayer. In this paper, we report 
the result from a study of the transverse sliding friction 
of a monolayer initially driven along the high symmetry 
direction of the periodic pinning potential. We find ev- 
idence of a transverse depinning transition at a nonzero 
temperature Tp. The nonequilibrium dynamic phase be- 
low Tp shows a critical depinning force in the transverse 
response even when the overlayer is in a state sliding 
smoothly along the longitudinal direction. 

We consider a model of interacting adatoms in a peri- 



odic potential subject to an external driving force [^,1). 
The dynamics is described by the Langevin equation 
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where rj is the adsorbate position, U ~ u{r.i) is the 
periodic substrate pinning potential, V = 
is the interacting potential between adatoms, F is the 
uniform external force acting on each adatom and is 
an uncorrelated stochastic force, with zero average, and 
variance related to the microscopic friction parameter rj, 
the mass of the particles m and the temperature T by 
the fluctuation dissipation relation < fi'{t)f"{t') >= 
2r]mkT5{t — t'), where a represents the vector compo- 
nents. We choose a periodic pinning potential with 
square symmetry 



u{r) — Uo[2 ~ cos{2-Kx/a) — cos(27ry/a)] 



(2) 



where a is the lattice constant of the substrate and an 
interacting Lennard- Jones pair potential 
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(3) 



where e is the well depth and rg is the particle separation 
at the minima in the pair potential. Periodic bound- 
ary conditions are used, with ro/a = 1.56 and the ad- 
sorbate coverage set to 9 = 1/2. With this choice, the 
ground state of the overlayer is a commensurate pinned 
c(2 X 2) structure which has been extensively studied in 
connection with the boundary lubrication problem 
and its longitudinal response to an applied force is well 
understood. Thus, this pinned structure provides a con- 
venient starting point for the present study of transverse 
response. We vary the temperature T and the ratio e/Uo 
that is essentially a measure of the stiffness of the over- 
layer relative to the pinning potential. Dimensionless 
units are used where a = 1, m = 1, Uo — ^ and we 
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normalize the force by 2TrUo and velocities by 21111 o/r]. 
We study the nonlinear response of the overlayer through 
dynamical simulations, using standard methods of Brow- 
nian molecular dynamics [T^ . As in the previous work 
on the longitudinal response small systems are used 
in order to allow for long equilibration times and reduce 
statistical errors. A system consisting of L x L substrate 
atoms with L = 10 to 20 was considered with the time 
variable discretized in units of 5t — 0.002 — O.OIt, where 
T = [ma? /UoY^"^ ■ Calculations were performed with typ- 
ically 1 — 4 X 10® time steps for each calculation of time- 
averaged quantities allowing equal time steps for equili- 
bration. 

As shown previously for an overlayer with the Lennard- 
Jones interacting potential 1^ and also for a model with 
pure elastic interacting potential |^ , the behavior of the 
drift velocity Vd as a function of F along the high sym- 
metry direction x or y depends strongly on the initial 
equilibrium phase at F = 0. We first summarize the 
known behavior for the longitudinal response along the 
high symmetry x-direction. When the overlayer at = 
is in an initial fluid or incommensurate state, then Vd will 
be nonzero for arbitrarily small external force leading to 
a finite sliding friction 1/fjL = Vd/F for arbitrarily small 
values of F. For our choice of parameters, the initial state 
has a commensurate c(2 x 2) structure below the melt- 
ing temperature T,„. For this initial commensurate state, 
the drift velocity is essentially zero below a critical value 
of F, neglecting the contribution from creep motion due 
to thermal activation. The F x Vd characteristics shows 
hysteresis behavior with unequal critical forces Fa and 
Fh for increasing or decreasing external force Fa,, corre- 
sponding to the static and kinetic friction forces, respec- 
tively. This behavior is shown in Fig. 1 for a temperature 
T = 0.2 much below the melting transition temperature 
Tm ~ 1.2 of the commensurate overlayer at = 0. In 
addition, at Fb there is a velocity gap Vf,. For sliding ve- 
locities V < Vb, smooth sliding is not possible. As a con- 
sequence, the motion of the overlayer shows stick and slip 
behavior when driven by a spring moving with a veloc- 
ity y < T4, in agreement with experimental observations 
where smooth sliding is only found above a critical slid- 
ing velocity [|l| . As first pointed out by Persson ||] , for a 
sufficiently strong interacting system, the particle veloc- 
ity distribution relative to the center of mass is Gaussian 
and one can define an effective temperature of the over- 
layer T* {kT* = m(< v'^ > - <v >^)/2). The irregular 
motion of the adatoms in the region Fb < F^ < Fc leads 
to T* > T, with the extra thermal dissipation, propor- 
tional to T* — T, being needed to balance the heating ef- 
fect of the external driving force to reach a steady state. 
The energy transfer to the overlayer is possible because 
the continuous translational symmetry is broken by the 
external pinning potential. In the limit F oo, the 
pinning potential is negligible, and the overlayer slides 
rigidly with T* — T. When the force is decreased, the 
sliding overlayer undergoes a dynamic melting transition 



into a sliding liquid phase at F^ where T* reaches a value 
above the melting transition temperature Tm- With fur- 
ther decrease of the external force, the overlayer is finally 
pinned at the critical kinetic frictional force Fb to become 
the commensurate phase again. 

In the sliding solid phase that exists at F^ > Fc in Fig. 
1, the motion of the overlayer essentially averages out the 
pinning potential. Therefore, for our choice of isotropic 
interacting Lennard- Jones potential, the overlayer should 
have a triangular structure although the pinning poten- 
tial has a square symmetry. However, in the transverse 
direction the pinning potential remains essentially static 
in a reference frame comoving with the overlayer [|l2| and 
it is plausible that the overlayer may remain pinned in 
the transverse direction, depending on temperature and 
velocity, although it is sliding in an incommensurate state 
along the longitudinal x-direction. 

We have determined the nonlinear response of the over- 
layer in the transverse direction in the sliding state. Once 
the system has reached a steady state for a fixed F^ > Fc, 
an additional force Fy is applied and the average drift 
velocity Vd is determined after the new steady state is 
reached. Alternatively, vt can also be obtained by apply- 
ing a total force {Fx,Fy) to an initial equilibrium state 
{Fx = 0). The two methods give the same result indi- 
cating a well defined transverse depinning behavior [T^ ] 
as shown in Fig. 2a. At temperatures below a critical 
value Tp, which is much smaller than the melting transi- 
tion temperature Tm, the mobility /it vanishes implying 
that the overlayer remains pinned along the transverse 
direction, even when it is sliding along the longitudinal 
x-direction. Above Tp, the mobility /it has a finite value. 
In the temperature range Tp < T < Tm , there is no trans- 
verse critical force and the overlayer is an incommensu- 
rate sliding solid in all directions. As shown in Fig. 2b 
and 2c, below Tp the Fy x Vy characteristics shows a crit- 
ical transverse force Fya and hysteresis behavior when 
the system is first equilibrated at lower temperatures. 
Note, however, that for decreasing transverse force, the 
pinned state is only reached at Fyb ~ 0, unlike the be- 
havior in the longitudinal direction in Fig. 1, suggesting 
that along the transverse direction the existence of a crit- 
ical force Fya may not necessarily lead to stick-and-slip 
motion. The behavior of the nonlinear response in the 
transverse direction correlates with the structure of the 
sliding state at Fy = 0. Fig. 3a, b and 3c, d show the 
instantaneous configurations and the trajectories of the 
adatoms in the overlayer below and above Tp for Fy — 0. 
The trajectories are obtained by superposing successive 
configurations in the sliding state for a fixed time inter- 
val. It can be seen from Fig. 3c and 3d that in the 
transverse pinned phase below Tp, the particles are es- 
sentially moving in one-dimensional channels along the 
principal x-axis of the pinning potential and the parti- 
cle coordination number is 4, leading to a square lattice 
structure. For T > Tp, the particle trajectories are free 
to fluctuate in the transverse direction (Fig. 3b) and an 
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isotropic triangular lattice structure occurs as shown in 
Fig. 3a. This is exactly the structure that the overlayer 
would have in the absence of a periodic pinning potential. 
The isotropic triangular structure has also been observed 
in previous calculations Q but was considered the only 
possible phase in the sliding state at any finite tempera- 
ture. Here we find instead that for temperatures below 
Tp a transverse pinned phase with different structure is 
the stable configuration. The dynamical phase transition 
separating the two sliding states has the characteristic 
features of an equilibrium first order transition. Besides 
the hysteresis effects mentioned earlier, we also observe 
supercooling effects in the sense that when the sliding 
state above Tp is cooled down to a temperature below 
Tp, at fixed F^, the final sliding state retains its trian- 
gular structure in a metastable state as shown in Fig. 
4a. This indicates that a nucleation process with large 
energy barrier is involved. While the details change, the 
hysteresis and supercooling effects are observed up to the 
largest size (20 x 20) studied here. We also note that the 
physics underlying the transverse depinning transition is 
very similar to that of the longitudinal depinning transi- 
tion which has been shown to be very general and model 
independent 

The results for other values of and e are qualita- 
tively the same as long as F^ > Fc- The linear transverse 
mobility nt ( in the limit i^y ^ ) as a function of the 
longitudinal driving force is shown in Fig. 4b for a tem- 
perature T — 0.2 below Tp. It is seen that /it remains 
zero for large F^ and the onset of finite /if for decreas- 
ing force occurs at F^ — Fc, where the sliding state goes 
from a solid phase into a fluid phase as would be expected 
from the behavior of the longitudinal response in Fig. 1. 
So there is a transverse depinning transition at all values 
of Fx > Fc- The dependence of the transverse depinning 
transition temperature Tp on the stiffness of the overlayer 
e together with the melting transition temperature Tm at 
Fx = are indicated in Fig. 4c. These transition temper- 
atures are located from the onset of finite linear mobility 
fit and fiL respectively for increasing temperatures. 

Recent work |To| ] on the dynamics of sliding two dimen- 
sional vortex lattices on a periodic potential also shows 
evidence of a transverse pinning and critical Lorentz force 
in the sliding state. The form of interaction and pinning 
potentials for the vortex lattice model are different from 
that for the overlayer model considered here, but this is 
not expected to be qualitatively important in the sliding 
state [Q. More importantly, the results for the vortex 
lattice were obtained at zero temperature and for an over- 
damped dynamics with zero vortex mass. Our results 
are obtained for a nonzero mass, finite damping coeffi- 
cient 77, and finite temperatures. We have shown that in 
our model, the transverse pinned state is stable against 
thermal fluctuations up to a transverse depinning critical 
temperature Tp < Tm as indicated in Fig. 4c. In view 
of the similarities of these systems and the current inter- 
est on the transverse critical current for vortex lattices 



in superconductors , it should be also interesting to 
investigate how a nonzero mass will affect the slid- 
ing state in vortex lattices and if the presence of weak 
disorder will broaden the sharp transition found in the 
present study. 

The existence of a transverse static friction in the slid- 
ing state implies that once two lubricated surfaces are 
sliding relative to each other in a particular direction 
commensurate with the underlying pinning potential, 
changing the direction of motion may require a nonzero 
additional force depending on the temperature and slid- 
ing velocity. Experimentally, the transverse critical force 
could in principle be detected from measurements of slid- 
ing friction with simultaneous transverse (in plane) vibra- 
tions between two boundary lubricated sliding surfaces, 
using similar techniques as in recent experiments which 
introduces coupling to normal (out-of-plane) vibrations 
and also in sliding colloidal layers confined between 
two surfaces Below the thermal depinning temper- 
ature the transverse vibrations should have very little 
effect on the longitudinal motion for sufficiently small 
vibration amplitude, while above the transition induced 
oscillations on the longitudinal friction are expected. We 
hope that this work will motivate novel experiments on 
sliding friction. 
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1.2- T = 0.2 




FIG. 1. Drift velocity Vx and effective temperature of the 

overlayer T* as a function of the external force along the 
a;-ax;is for e = 1 and rj = 1. The direction of variation of Fx 
is indicated by arrows. 



FIG. 2. (a) Transverse mobility as a function of temper- 
ature obtained by applying a total force {Fx,Fy) to an initial 

equilibrium state (squares) and by applying and additional Fy 
to the sliding state at {Fx,0) (circles), (b) and (c) Transverse 
velocity Vy and effective temperature of the overlayer T* as a 
function of an additional force along the j/-axis for fixed Fx- 
For large Fy, the behavior of T* is smilar to Fig. 1. Large 
arrows indicate the direction of force variation. Results are 
for L = 10, Fx = 1.5 and e = 1. 




FIG. 3. Snapshot pictures ((a) and (b) ) and trajectories 
of the adsorbates ((b) and (d)) at Fj, = for temperatures 
T = 0.4 and T = 0.2 above and below the depinning tem- 
perature Tp, respectively, for Fx — 1.5 and e = 1. Grid lines 
represent the periodic pinning potential. 
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(a) T=0.2 
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FIG. 4. (a) Snapshot picture of the final state at T = 0.2, 
below Tp, obtained from an initial sliding state at T — 0.4, 
above Tp, by lowering the temperature at fixed Fx = 1.5. (b) 
Transverse mobility fit as a function of the longitudinal force 
Fx at a temperature below the depinning transition temper- 
ature Tp. Results are for e = 1, L = 10. (c) Critical tem- 
perature Tp for transverse depinning in the sliding state for 
Fx = 1.5 and equilibrium {Fx = 0) melting transition tem- 
perature Tm as a function of the overlayer stiffness e. 
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